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In this note we study the linear dynamics of scalar graviton in a de Sitter background in the 
infrared limit of the healthy extension of Hof ava-Lifshitz gravity with the dynamical critical exponent 
z — 3. Both our analytical and numerical results show that the non-zero Fourier modes of scalar 
graviton oscillate with an exponentially damping amplitude on the sub-horizon scale, while on 
the super-horizon scale, the phases are frozen and they approach to some asymptotic values. In 
addition, as the case of the non-zero modes on super-horizon scale, the zero mode also initially 
decays exponentially and then approaches to an asymptotic constant value. 
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I. INTRODUCTION 

Recently, a power-counting renormalizable ultraviolet (UV) complete quantum gravity theory was proposed by 
Hofava [1, 2]. This theory is characterized by the anisotropic scaling between time and space, so the complete 
diffeomorphism invariance of general relativity (GR) is lost, instead the Hofava-Lifshitz (HL) gravity is invariant 
under the so-called "foliation-preserving" diffeomorphism Diff(Af, F). Since this theory was proposed, a great deal 
of efforts have been made, including studies of cosmology [3-17] and black hole physics [18-25], among others [26- 
28]. Because of the differences of diffeomorphism groups between in HL and in GR, one expects to see some new 
dynamical degrees of freedom of gravitational fields in HL gravity. In the minimal realization [2] , because the lapse 
function N(t) introduced in the Arnowitt-Deser-Misner (ADM) formalism respects the "projectability condition", 
i.e. it is a function of time only, the usual local Hamiltonian constraint becomes into a global one which will not 
affect the local dynamics. So, the absence of local Hamiltonian constraint leads to a new scalar degree of freedom 
in addition to the usual helicity-2 polarizations of the graviton at the linear perturbation level [2, 29-38] (see [39] 
for a review on the "projectable model"). However, further studies on the non-linear dynamics show that the extra 
scalar mode suffers several pathologies, such as instability and strong coupling problem [40, 41]. A naive extension 
to the Hofava's original proposal, which is often refereed to the "nonprojectable model" in literatures [1, 2], is to 
restore the full dependence of the lapse function on space and time N(t, x). In this model the scalar graviton becomes 
non-dynamical, because the equation of motion for the lapse function gives the local Hamiltonian constraint [42, 43]. 
Unfortunately, the "nonprojectable model" also confronts several conceptual difficulties [41, 44-46]. Recently, a new 
"projectable extension" with gauged (7(1) symmetry is proposed by Hofava and Melby- Thompson [47]. 

In this note we will investigate another extension of the HL theory, which is called healthy extension of HL 
gravity in literatures [41, 48] (see a nice review in [49]). In this model, a new ingredient = VjiV/JY, which 
transforms under Diff(Ajf, T) as a spatial vector and a time scalar, is introduced into the action and the Hamiltonian 
constraint becomes into the second-class one. As a result, one extra degree of freedom should appear in the healthy 
extension. A preliminary parameterized post-Newtonian study in [48] shows that the extra scalar graviton is free 
from ghost instabilities in some parameter regions. Furthermore, the cosmological evolutions of metric and density 
scalar perturbations in both radiation and matter dominated eras are investigated in [50]. They found that, although 
the system has two scalar degrees of freedom, corresponding to a scalar graviton and an adiabatic matter fluctuation, 
the late-time evolution of perturbations can be sufficiently specified by the value of one gauge-invariant variable. So, 
it is natural to ask how the scalar graviton behaves in the early universe, i.e. de Sitter or inflationary phase. In 
this note we study the linear dynamics of HL scalar graviton in a flat universe with a positive cosmological constant. 
Both our analytical and numerical results show that the non-zero Fourier modes of scalar graviton oscillate with an 
exponentially damping amplitude on the sub-horizon scale, while on the super-horizon scale, the phases are frozen 
and amplitudes continuously decays until they approach to their asymptotic constant values. In addition, as the case 
of the non-zero modes on super-horizon scale, the zero mode also decays exponentially with respect to time initially, 
and then approaches to an asymptotic constant value. 
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This note is organized as follows. In section II, we briefly review the "foliation-preserving" gauge symmetry and the 
setup of the healthy extension of Hof ava-Lifshitz gravity in 3 + 1 dimensions. Then we construct some gauge-invariant 
variables and derive the background equations in a spatially fiat universe in section III. In section IV, we present 
both analytical and numerical studies on the cosmological linear perturbations. Finally, we conclude in section V. 



II. THE HEALTHY EXTENSION OF HORAVA GRAVITY 



In this section, we firstly present the gauge transformations compatible with the "foliation-preserving" diffcomor- 
phism Diff(M, J 7 ), then briefly review the setup of the healthy extension of HL gravity in 3 + 1 dimensions. 

A. Gauge symmetry 

The field contents in the healthy extension of HL gravity are 

lapse: N(t,x);, shift: 7V*(t, x) , 3d spatial metric : gij(t,x) , (1) 

where we abandon the "projectablility condition" on lapse function N(t,x). In terms of these fields we can cast the 
4-dimensional line element in the ADM formalism as 

ds 2 = -N 2 c 2 dt 2 + gij [dx l + N'dt^j (dx j + N j dt) , (2) 

where x° = ct and the light speed c is restored explicitly in order to obtain the nonrelativistic theory from relativistic 
one by taking c — > oo limit. 

In the Hof ava's proposal the local Lorentz invariance is violated due to the anisotropic scaling between space and 
time. For instance, in 3 + 1 dimensions the coordinates (t, x) scale as 

t->t z t , X4<x, (3) 

where z is called dynamical critical exponent. In the case of general z, the classical scaling dimensions of the fields 
are 

[N i ] = [c]= [ ^ = z-1, [ gij ] = [N}=0. (4) 

In order to investigate the spacetime diffcomorphisms in the nonrelativistic theory, it is convenient to start with 
the relativistic metric g^ v but with c restored 

/ -N 2 + N i N i /c 2 , Njc\ 
9tw ~ { Ni/c, 9ij ) ' (5) 

( - 1/N " ' N*N~ 2 /c \ 

9 ~ \ N l N- 2 /c , gv - N l NiN- 2 /c 2 ) ' ( ' 

where we only keep the leading terms in 1/c expansion. Correspondingly, the covariant generators a: ^ = + ^ of 
spacetime diffeomorphisms are also expanded formally with respect to the small parameter 1 /c 

e = cf(t)+o(i/ C ), e = c(i,x)+o(i/c 2 ). (?) 

As stated above, the nonrelativistic transformation rules are easily obtained by taking c^oo limit of the relativistic 
diffcomorphisms 

5gij = - 5jfeViC fe - gtk^ 3 C k - fkj , (8) 

SNi = -VitfNj-eVjNi-pgij-fNi-fNi, (9) 
SN = -( j VjN - fN - fN , (10) 

where Vj is the 3-dimcnsional covariant derivative compatible with spatial metric g^ and the over dot denotes the 
derivative with respect to the nonrelativistic time t. 
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B. Model setup 



Comparing with the naive promoting the lapse function to a spacetime field in the "nonprojectable model" [1, 2], 
a new ingredient for constructing gauge-invariant terms in the action is introduced in the healthy extension of HL 
gravity [41, 48] 

a t = VilniV(i,x) . (11) 

In fact, the appearance of this term in the action is compulsory to make theory be free of the pathologies, which has 
been mentioned in the previous section. The key point of the healthy extension is that once terms with this new 
ingredient are introduced in the action, the Hamiltonian is no longer linear in the lapse, and the equation of motion 
for N(t, x) becomes into the second-class constraint. By using the standard Hamiltonian analysis [51, 52], the number 
of degrees of freedom in the D + 1 dimensional healthy HL gravity is 

JV=^(dim7>-2Ci-C 2 ) = ^D(D-l), (12) 

where dim V = (D + 1)(D + 2) is the dimension of phase space, C\ = 2D is the number of first-class constraints, and 
C2 = 2 is the number of second-class constraints. So, in the 3 + 1 spacetime this model exhibits an extra degrees of 
freedom (dof) in addition to the usual transverse traceless helicity-2 gravitons. 

Because our primary purpose is to investigate the linear dynamics of scalar perturbations on a flat universe, we 
write the most general form of the action relevant to our interests 



J dtd 3 x^N{c K - Vfotf.Oi]} , (13) 



where the kinetic term takes usual form 



Ik = KijK» - XK 2 , (14) 



with the extrinsic curvature defined as Kij — (gij — ViNj — \7 jNi)/2N . The most general potential up to the terms 
with dynamical critical exponent z — 3 is given by V\g%j, a,] = Y^ z =i ^ z > 

Vi = —g\R — riciiO 1 + a , (15) 
V 2 = M p - l 2 (g 2 R 2 + .7.., /.',,/.•" + m^Aa* + rfyRV t a % ) , (16) 
^3 = M^faRAR + grXJljiVH' 1 + r^AV + ^ARV.a 1 ) , (17) 

where Rij and R are the Ricci tensor and Ricci scalar, respectively. In our convention the Laplace operator reads 
A = ViV 1 , its square and cubic are A 2 = A • A and A 3 = A • A • A. A parity-violating term e^ k RuW jR l k is absent 
in our paper because it only affects the tensor cosmological perturbations at the linear order. However, this term 
will contribute to the scalar perturbations at the high-order level due to the appearance of coupling vertices between 
scalar and tensor modes. 

Now we would like to write out the equation of motion for this system of gravity. Variation with respect to the 
lapse N(t,x) gives the local Hamiltonian constraint 

C K + V + 277VX - ll-AVX + -J* AR ^A 2 V lfl - + ^A 2 R = . (18) 
lvl P i m P i m P i m P i 

Equation of motion (eom) for Ni gives the ordinary momentum constraint 

VjTT ij = , (19) 

with 7r y = K.3 % — \Kg % 3 . In principle, the propagating equations of this system can be obtained by the variation of 
gravitational action with respect to g^, i.e. SSu/Sgij = 0. However, the explicit expressions of eom for g^ are very 
lengthy, and one can find them in [50] . In the section IV, we will derive the linear approximations of these equations 
through a perturbative approach. 
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III. GAUGE-INVARIANT VARIABLES AND BACKGROUND EVOLUTIONS IN A FLAT UNIVERSE 

In this section we will firstly study the "foliation-preserving" gauge transformations of linear scalar modes in a flat 
universe, and then demonstrate that the dynamics of scalar perturbations at the linear level can be analyzed without 
complete gauge fixing due to the novel properties of "foliation-preserving" gauge transformations, finally derive the 
set of background equations by using the linearized action. 

The background line element of a flat Friedmann-Robertson- Walker (FRW) universe reads 

ds 2 = -dt 2 + a 2 5 ij dx l dx j , (20) 

where the scale factor a(t) is a function of time only. In our convention we linearize the metric scalar perturbations 
as 

$9oo = -20(i,x) , 
Sgoi = a 2 diB , 

5 gij = -2a 2 (^<% - didjE) , (21) 

where di is the ordinary spatial derivative compatible with Kronecker delta function <5 y and d 2 = d l di. The transfor- 
mation rules under "foliation-preserving" diffcomorphism for these modes can be obtained by virtue of (8)-(10), 

54> = -f{t) , SB = -e(t, x) 5i> = f(t)H(t) , SE = -e(t, x) , (22) 

with £ J = die + e 1 and die' 1 = 0. 

Armed with these transformation rules, we construct some useful gauge invariant variables as follows 

* = 0+(|)\ (23) 
= B-E, (24) 

* = ip + H [ <f)dt, (25) 



J t 

as pointed out in [30], ^ in fact is the integral version of $. 

From (22), we can see that the gauge variations for the quantities <j>(t,x) and tp(t,x) are two functions of time only. 
This novel feature of gauge transformations allows us to analyze the dynamics of linear perturbations by using the 
gauge dependent variables, which will not lead to any gauge artifacts even though we do not fix the gauge completely. 
To see this more clearly, here we give a simple example to illustrate our method. We consider a system consisting of 
two fields x an d 7 

s\xn\ = -\ J {x 2 - dixd'x- m 2 x 2 - c [c) [xn\} , (26) 

where we explicitly write down the free sector for x, but denote both the 7 sector and interaction terms by an implicit 
form £( c \ Furthermore, we assume that the action (26) is invariant under the transformation 

X(i,x) -> x(i,x) =x(*,x)+A(i) , (27) 
7(t,x) -> 7(t,x) =7(t,x)+fl(i) . (28) 

Under the "gauge" (^4 = 0, B = 0), the eom for x reads 

X-d 2 x + m 2 x + C[x,7] =°> ( 29 ) 

with the coupling term C = ^SC^/Sx- On the other hand, we can also write the action (26) in terms of tilde 
quantities in (27) and (28) as 

S[ X , 7 ] = S[x,7] = ~l J {x 2 -^x-m 2 x 2 -£ (c) [x,7]} , (30) 
= S[ X , 7] - \ J 2(ix - m 2 A X ) -\J [X, l,A,B]-\ J F[A, B] , (31) 
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with 

[x, 7] = ^ (c) \X, 7] + [x, 7, A, B] + ^ (2) [4 B] , (32) 
T[A,B] D C {2) [A,B] , (33) 

where [x, 7, -A, -B] contains the linear terms of \ om y, while [A, B] and J"L4, B] contains the quadratic terms 
of A, B only. Thus, the variation of tilde action (30) with respect to \ gives 

X-d 2 x + m 2 x + C{ X ,~f}=S[A,B} , (34) 

with S = —^SjC^/5x — A — m 2 A. Comparing with the result in (29), equation (34) has an extra source term on 
the right hand side, but it depends on time only. Given this fact, if we turn to the Fourier space, we can easily see 
that such only time dependent source term vanishes for all k 7^ modes, while it contributes a Dirac delta function 
to k — mode. So, we conclude that the only time dependent gauge transformations do not affect the dynamics of 
k 7^ mode, i.e. the eom for k 7^ modes is the same in all gauges. These features of Hofava-type theory allow us 
to study the dynamics of non-zero Fourier modes by using the gauge dependent variables </> and ip directly without 
leading to any gauge artifacts. 

In the rest part of this section, we list the background equations which can be obtained by varying the linearized 
action with respect to the scalar perturbations. In details, eom for <j> gives Hamiltonian constraint 

a + 3(1 - 3X)H 2 = , (35) 

eom for ip or d 2 E gives the evolution equation 

H = 0. (36) 

Finally, the momentum constraint, i.e. eom for diB, is trivially satisfied on the background. Note that (35) and (36) 
give a de Sitter solution. 



IV. DYNAMICS OF SCALAR MODES IN LINEAR COSMOLOGICAL PERTURBATIONS 



In this section we will investigate the linear dynamics of scalar perturbations without any gauge fixing. First, we 
have to derive the second order action 



5 (2) = K 
H 2 



dtd 



>x a 3 j(l - 3A) 3(V> + H<j)f + 2(tp + H<j))d 2 {B - E) 



+ (1 - X)d 2 (B - E) ■ d 2 (B -E)- ya' 2 (j)d 2 (/) + 2g ia - 2 (24> - *p)d 2 ip 



(16ff 2 + 6# 3 )d 2 V> • d 2 i> - m d 2 4> ■ d 2 (f> + 4rj 3 d 2 ip ■ d 2 
(I604 - 60 5 )<9 2 V • <9 4 V> - md 2 4> ■ d A (f> + 4r/ 5 <9 2 • d 4 V>] } 



(37) 



where we have used the background equations and dropped all surface terms. Hamiltonian constraint is given by eom 
for <j>, i.e. S^sff = 



(1 - 3A) 6H 2 <j> + 2Hd 2 (B - E) + 6Hip - 2f] a - 2 d 2 <p + Ag^d 2 ^ 
+2M p| 2 ))2 a- 4 9> - AM p 2 m a- A d^ + 2M pl 4 ma - 6 d 6 (b - AM^a^d^ = . 

The eom for d 2 (B — E) produces the momentum constraint 

(1 - 3\)(H<j> + ip) + (1 - X)d 2 (B -E) = 0. 



(38) 



(39) 



Both the Hamiltonian and momentum constraint equations are consistent the linear approximations of (18) and (19). 
Finally, the propagating equations is obtained by eom for tjj 



(1 - 3A) - 6Hd 2 (B -E)- 2d 2 (B - E) - 18H*<f> - 6H<j) - 6H<f> - 18i?V - 6ip + Agia^d 2 ^ - i>) 



2(16 ff2 + 6 53 )5 4 V' + 4r ? 35 4 



2(16 ff4 - 6 95 )5V + 4%3 6 



. 



(40) 
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A. Non-propagation of the scalar graviton in the original non-projectable version 

In order to check the method presented in the section III, we need to recover the results in the original "nonpro- 
jectable model" by using this method. So, we turn off all a 1 terms in action (13) and (37). For further simplification, 
we turn off all operators with dimension higher than 2 in the gravity potential (16) and (17), i.e. merely keep R and 
a terms in the potential (15), without losing any information. Thus, the second order gravity action (37) becomes 

M 2 fry 

S h } = -yW dtd 3 x a 3 {(l-3\) 3(ip + H(j)) 2 + 2(ip + H(j))d 2 (B - E) 

+(1 - \)d 2 {B - E) ■ d 2 (B -E) + 2g ia - 2 {2<t> - V)<9 2 v} ■ ( 41 ) 
Solving the Hamiltonian and momentum constraints gives 

♦--B^rOT**- (42) 

ffii = 9i^-d 2 f , (43) 
and plugging them back into (41), we arrive at 

S%> = ^Lj dtd 3 x a 3 ^ gia - 2 p 2 ^ - 2 9l a' 2 ^ 2g 2 ^±^d^ ■ d 2 ^ . (44) 

From the above expression we can see that the kinetic term ip 2 is absent, so at the linear perturbation level there does 
not exist the dynamical scalar graviton in the "nonprojectable model" . This result agrees with the observations made 
in [42, 43]. Beyond the non- linear order, however, the conclusion is changed. As shown in [47], this mode becomes 
a dynamical one if the terms of the fluctuation beyond the quadratic order are taken into account in the action. In 
addition, the absence of the usual quadratic time derivative term leads to an unacceptably strong-coupling problems 
on the small scale [53] . 



B. Linear dynamics of scalar graviton in a de Sitter universe 



As be expected from Hamiltonian analysis (12), there should exist one extra scalar graviton in the pure gravitational 
system for the healthy extension of the HL theory. Now, we will study the linear dynamics of scalar graviton in a flat 
de Sitter universe with the healthy extension of the HL gravity. 

Solving the momentum and Hamiltonian constraints formally, we obtain 



d 2 (3 



1-3A 
1 - A 



(iP + Heft , 



and 



2 (-^d t - 9l a- 2 d 2 + M- 2 V3 a- 4 d 4 + AC^o" 6 ^ 



$ - r)a- 2 d 2 + M'^a^d 4 + M^^a^d^ 



(45) 



(46) 



Firstly, we plug the momentum constraint into action (37) and reduce it into a functional of two scalar fields <p an d 4> 

2(1 - 3A) 



S 



f = ^kf dt s x fl 3 



1-A 



-{ip + H4>Y - V a 



-2 j.a2. 



+ 2g ia - z {2<j)-ij)d 2 i> 



-M~ 2 a-* 



(16 5 2 + 6 53 )5 2 V- • d 2 ^ - r, 2 d 2 4> ■ d 2 <j> + 4773d 2 V- • 3 2 <; 
(I654 - 6g 5 )d 2 ip ■ <9 4 V - md 2 4> ■ d A (t> + Ar] 5 d 2 4> • d 4 V } 
where <j) ne ld is non-dynamical dof and its eom provides a constraint 

I _ Q\ 

m{^ + H4>)-2f 1 a- 2 d 2 4> + Ag 1 a- 2 d 2 iP + 2M- 2 m a- 4 d 4 q 



(47) 



1-A 



-4M o - 2 77 3 a- 4 dV + 2M p - 4 r/ 4 a- 6 a 6 - AM^a^d 6 ^ = 



(48) 
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The eom for ip is 

4(1 -3A) 



1 - A 

V 



$ + 3H(ip + H4>) + H<P + H(f> + 4 9l a- 2 d 



V) - M , 2 a- 4 2(16 ff2 + 6 53 )5V + 4r ?3 3 4 < 



-M-, 4 a~ 6 



2(16 ff 4-6. 95 )3 ( V + 4?? 5 c> 6 



. 



(49) 



In the Fourier space, the above two partial differential equations reduce into the ordinary one 

1 - 3A 



1-A 



■AH(ip + H<P) + 2r]a- 2 k 2 (f> - Ag^k 2 ^ + 2M- 2 r] 2 a- 4 k 4 (f> 



and 



-4M- 2 m a' 4 k 4 ^ ~ 2M- l 4 m a- 6 k 6 <j ) + AM^a' 6 ^ = 
4(1 -3A) 



(50) 



1-A 

-M~ 2 a- 4 k 4 



ip + 3H(ip + H<f>) +H(j) + H<j) - A gi a- 2 k 2 (<f> - i/>) 



2(16 ff2 + 6.g 3 )V' + 4?73(i 



2(16g 4 - 655)7/" + 4?75<; 







(51) 



where we have suppressed the momentum k implicitly in the Fourier modes 4>k and ipk- Before solving the above 
differential equations, we would like to discuss the ghost-free condition which comes from the positiveness of kinetic 
term ip 2 in the Lagrangian (47). Because the field </> is not a dynamical one, in order to analyze the dynamics of the 
extra scalar graviton easily we need to solve the Hamiltonian constraint equation (50) and plug it into the Fourier 
form of the action (47), then the action becomes a functional of tp only. In the following discussions on the ghost- free 
condition we will take this approach. Firstly, let us focus on the ultraviolet limit (fc — > 00), in which fc 4 ' 6 terms are far 
more important than k ' 2 terms. Because there are some coupling constants (772, %, f]i, %) of higher spatial derivative 
terms, he ghost- free condition can be easily satisfied, i.e., the parameter constraints to avoid ghost instability are 
rather loose, here we will not write them down explicitly. On the other hand, the ghost-free condition in the infrared 
limit (k — > 0) is controlled by the parameter 77 only, one can obtain an explicit constraint on 77. In the infrared limit, 
we can safely neglect 0(k 4 ' 6 ) terms, thus (50) reduces into 



4^1 k 2 ip - ^Y^AHtP 
^4H 2 + 2r)k 2 



1-A 



H 



(l-A)gifcV 
(1 - 3A)tf 2 2(1 



(1 - \)r)k 2 ip 
3A)i? 3 



(52) 



where we have denoted the physical wavelenght by k — k/a. The relevant part for the kinetic term in the action reads 



? (2) 



M 2 pl 
2(2tt) 3 



dt / d 3 ka 3 



2(1 -3A) 
1-A 



\ijj + H<b\ 2 + k 2 V \<f>\ 2 



(53) 



Here we emphasize that the Fourier components <f>, ip are complex fields, in order to keep the action real we have the 
absolute value symbol for each field component, or we can also write \(j>\ 2 as <f>{t, k)</>*(i, k). Plugging (52) into (53), 
one can easily see that the leading term in the kinetic term reads 



Kinetic term : 



r]k 2 



(54) 



from which we can see that the ghost-free condition in the infrared limit is 77 > 0. Our result agrees with the 
stability constraint in the Minkowskian background [48], but contradicts with the one in [55] where by Cerioni and 
Brandenberger find the ghost-free condition is 77 < 0. At the end of the paper, we will have more to say on this point. 

Due to the cosmological interests, in the rest of this section we will investigate both analytically and numerically 
the dynamics of scalar graviton in the infrared limit, i.e. we turn off all coupling constants (772, <?2: ^37 93, • • •) i n the 
front of the operators with dimension larger than 2. Firstly, we will solve the above equations numerically with the 
initial condition ip(to,'k) = 1, dip/dN\ to M — 0, (see Fig. 1). In the numerical calculations, we plot the dynamical 
evolutions of six different Fourier modes which cross horizon (fc» ~ a*H) at the efolding number = ln(a*/ao) = 1 
(Top left), TV, = 2 (Top right), N r = 3 (Middle left), = 4 (Middle right), TV* = 5 (Bottom left), and N* = 6 
(Bottom right), where star symbol represents for the moment of crossing horizon and the subscript zero for the initial 
time of de Sitter phase. The numerical results show that the Fourier modes oscillate with damping amplitudes on the 
sub-horizon scale, while freeze to some asymptotic values on the super-horizon scale. 
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FIG. 1: This figure shows the infrared behaviors of different Fourier modes of the scalar graviton if) in a de Sitter universe, 
where the horizontal axis denotes the efolding number TV = ln(a/ao). The six curves are the Fourier mode which exits Hubble 
horizon (fc* ~ a*H) at the efolding number JV. = 1 (Top left), N, = 2 (Top right), N, = 3 (Middle left), N, = 4 (Middle 
right), JV* = 5 (Bottom left), and N* = 6 (Bottom right), respectively. In our numerical calculations, we set the initial scale 
factor ao = 1, Planck mass M v i = 1, cosmological constant a = 1, gravitational coupling parameters A = 1.05, r\ = 0.1 and 
<7i = 1. The initial conditions are given by ip(to, k) = 1 and dip/dN\t 0> k = 0. 



Secondly, we will try to find some analytical solutions which are able to explain the numerical behaviors. For this, it 
is convenient to use the conformal time r = — (aH)^ 1 . Combining (50) with (51), we obtain the propagating equation 
for the scalar graviton in conformal time 



Fx (Hkr) i>" + F 2 (Hkr) V + F 3 (Hkr) \ifi = , 



(55) 



where prime denotes the derivative with respect to conformal time ' = d/dr. Notice that the equation (55) is 
obtained non-perturbatively, and the expressions for coefficients F\ F 2 and Fs are given in Appendix A. From the 
above expressions we can see that, the coefficients F\ , F 2 vanish if all at terms are absent in the action, i.e. rji = 0. 
So the solution of (55) is trivially ip = 0, which means that in this case the scalar graviton does not exist in the case 
of absence of matter sources. This result is consistent with our observations for the "nonprojectable model" in the 
subsection IV A. 

Now, we will turn to the healthy extension model in which not all coupling constants rji vanish simultaneously. 
Because the dimensionless quantity fcr is much less than unit on the super-horizon region, we can use it to expand the 
equation (55) pcrturbatively, and solve them order by order. For the super-horizon modes {kr <^ 1), the propagating 
equation (55) reduces to 



tp" + [ik 2 Ti\>' + c 2 s k 2 i/> = 



with 



M = 



1 - A 
3A- 1 



V ■ 



(56) 



(57) 
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and sound speed 



A - 1 
3A- 1 



v 



(58) 



where we have neglected the contributions from coupling terms r]2 , 52 , V3 1 93 etc The dynamical stability can be 
satisfied provided the sound speed is positive definitely (c 2 > 0) . Considering the constraint on the coupling constant 
(0 < r] < 2) from Minkowskian space [48], such stability condition can be achieved in the parameter regions A e 
(-oo, 1/3) U(l, +00) and 9l e (-00, -1) |J(0, +00). 

In order to compare with our numerical results, we rewrite the equation (56) in terms of the efolding number 
N = ln(a/a ) 



d 2 ip 
dN 2 



i H 2 



dip 
dN 



cjk 2 
a 2 H 2 



<P = o, 



(59) 



where the initial scale factor cto can be set to unit. For super-horizon modes, k/aH <C 1, the above equation reduces 
into a simple form 



d 2 tp 
dN 2 



<hp_ 
dN 



. 



The general solution of the above equation reads 



V> = de~ N + Q 



(60) 



(61) 



with two integrate constants C\ and C2 . From the above solution, we can see that the phases for super-horizon modes 
are frozen and the amplitudes experience an exponential suppression initially, and then approach to some asymptotic 
constant values. These behaviors are in agreement with our numerical calculations. 
For the sub-horizon mode (fcr » 1), (55) becomes into 

2 



tp" - -ip' + k 2 c 2 s iP = , 



with another sound speed 



A - 1 
3A- 1 



fh 



(62) 



(63) 



To avoid the exponential instability, the parameters must satisfy A e (—00, 1/3) (J(l, +00) and g\ e (—00, 0) (J(l, +00). 
The solution for (62) reads 



r 3/2 



C\ cos(c s fcT) 



sin(c s A:7 
c s kr 



C 2 I sin(c s fcr) + 



cos(c s fcr) 
c s kr 



where C\, C2 are another two integration constants. If we rewrite the equation (62) with the efolding number 

d 2 tp dtp c 2 k 



dN 2 



3 ^ + ^ = °- 



Then the general solution can be expressed as 



c s k 



■ cos 



c s k 
aH 



c s k 
~a~H 



+ C 2 



cos 



c s k 
~a~H 



c s k , 



c s k 
~a~H 



(64) 



(65) 



(66) 



From (66) we can easily see that the sub-horizon Fourier modes have two oscillating solutions. The amplitude of one 
solution is large, but it is decaying with the scale factor a -1 ; the amplitude of the other is constant, but it is smaller 
than the former by a factor aH/c s k <C 1. The Middle right, Bottom left and Bottom right panels in Fig. (1) show 
that the former exponentially suppressed modes dominate the sub-horizon solutions in our numerical calculations. 

Finally we would like to make a comment on the zero mode, for which the above analysis becomes invalid. Fortu- 
nately, we can rewrite the action (47) into a gauge-invariant form, if we neglect all spatial gradient terms by using 
the gauge-invariant variable ^ (25) 



c( 2 ) ~ 

a H — 



2 

M 2 pl 



dtd 3 x a 3 
dtd 3 x a 3 



2(1 - 3A) 
1-A 



(V- + H<j>f 



2(l-3A) ip 



(67) 
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Consequently, eom for ^ is given by 

* + = , (68) 

with the general solution 

* = Ae- 3Ht + B . (69) 

we see that as the case of the super-horizon non-zero modes, the general solution of also contains two parts, one of 
them decays exponentially with respect to time, the other one is a constant. 



V. CONCLUSION AND DISCUSSION 



In this note we investigated both analytically and numerically the linear dynamics of scalar graviton in the healthy 
extension of the HL theory in a de Sitter background. We found that due to the "foliation-preserving" diffeomor- 
phism Diff (M, F) , the gauge transformations of metric perturbations <j> and tp are two functions of time only. These 
novel features of gauge transformations allow us to analyze the dynamics of linear perturbations by using the gauge 
dependent variables, which will not lead to any gauge artifacts even though we do not fix the gauge completely. Given 
these observations, we studied the linear dynamics of scalar perturbations without any gauge fixing in section IV. 
Firstly, we used our method to study the dynamics of scalar graviton in the original "nonprojectable model" and 
found that the scalar mode does not exist in the absence of matter sources, which is consistent with the existing results 
in the literatures. Secondly, we analyzed the dynamics of non-zero Fourier modes of scalar graviton in the infrared 
limit of the healthy extension of HL gravity by our method. Both our numerical and analytical solutions show that 
on the sub-horizon scale, the non-zero Fourier modes oscillate with exponentially damping amplitudes; but on the 
super-horizon scale, the phases are frozen and amplitudes decay continuously until they approach to their asymptotic 
values, i.e. the Fourier modes are conserved once they cross the horizon. Finally, the dynamics of zero Fourier mode 
of scalar graviton is also presented by using the gauge- invariant variables. As the case of the non-zero modes on 
super-horizon scale, the zero mode also decays exponentially with respect to time initially, and then approaches to an 
asymptotic value. 

In this note we have only investigated the dynamics of scalar graviton in a de Sitter universe in the healthy extension 
of the HL theory without matter sources. For a more realistic case, we should take some matter sectors into account. 
In particular, in the early universe, it is natural to introduce other scalar field(s) to drive the inflation. Thus, in such a 
system there will exist two scalar perturbations at least, one is the HL scalar graviton and the others are the fluctua- 
tions of matter fields. So, it is natural to ask among all scalar perturbations which one is responsible for the generation 
of primordial adiabatic perturbations. Some other cosmological aspects, such as primordial non-Gaussianities and 
trans-planckian problems, in the healthy extension of HL gravity are also worthy to further investigate. 

Note added: At the last stage of this work, two papers [54, 55] appeared in the arXiv, which contain some 
relevant discussions. The authors of both [54] and [55] consider a system within the healthy extension of HL gravity 
coupled minimally with a Lifshitz-like inflaton field tp. In such a system, there exist two dynamical scalar degrees of 
freedom on the linear perturbation level around a cosmological background, one is the scalar graviton ip which newly 
appears in the healthy extension, the other is the inflaton fluctuation Sep. The authors in both [54] and [55] found 
that the extra scalar graviton couples with the matter field in a general FRW universe. But in the infrared limit, 
these two modes decouple. However, the reasons given in the two papers for the decouple phenomenon is a little 
bit different. In [54] Koh and Shin argued that the prefactor of kinetic term ip vanishes at the leading order 0(k ) 
when the momentum (k — > 0), i.e. the extra scalar graviton becomes non-dynamical due to the absence of kinetic 
term. On the other hand, Cerioni and Brandcnbcrgcr [55] performed the calculation to the next order and they found 
that the kinetic term reads k 2 ip 2 /H 2 at 0{k 2 ) order, but the mass of the extra scalar graviton becomes infinitely 
heavy when one recast the kinetic term into the canonical form (ip 2 /2), because the mass term takes the form of 
(mofc + m 2 fc 2 + mJz A + ■ ■ -)H 2 ip 2 jk 2 . Although the kinetic term does not vanish at the sub-leading order, the scalar 
graviton still becomes non-dynamical due to the infinitely heavy mass. Hence, at late times, it will decouple from low 
energy physics and will not contribute to cosmological perturbations on scales relevant to current observations. In 
our present work, we merely considered a pure gravity theory and found that the scalar graviton is still dynamical in 
the infrared limit. This result is consistent with those in [54, 55], because when one turn off the matter sector related 
with the inflaton ip in the Lagrangian, both mo and m-i vanish simultaneously, and the leading term becomes m4k A ip 2 . 
Furthermore, the kinetic term still takes the form of k 2 ip 2 /H 2 , so after expressing the Lagrangian in the canonical 
form the mass term of ip reads mJi 2 H 2 ip 2 , i.e. the extra scalar graviton keeps massless in the infrared limit of the 
de Sitter background. This feature is the most important difference from the theory coupled with scalar matter fields 
[54, 55]. Finally we would like to mention that the ghost free condition in a de Sitter background in our derivation 
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(rj > 0) are in contradiction from those obtained by Cerioni and Brandenberger (77 < 0) [55]. However, our result 
is consistent with the one from the consideration in a Minkowskian background [48]. So, it is necessary to fix this 
disagreement, but it is out of the scope of the present work. 
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Appendix A: 



In this appendix we list the coefficients in (55), which is obtained through a non-perturbative method. As a general 
consideration, we keep all terms in the action (47) 



2aM 2 pl 



F 2 (x) = 



8<j 
9(1-A) 



rG- 1 (x) 



3(1- X)J ' 

- VAX 6 + (4*MZ - 3 (1 - A)) V2 x 4 + (2aM p l l - 3 (1 - A) M 2 t ) r,x 2 



F 3 (x) = 8G- 1 (x) M-frjlx 12 - 2M pl %r l5 x w + M p 2 {r, 2 - 2g lV6 ) x s + (^2g im + M p f ^ 



V5 x 



+8a (1 - Xy 1 G- 2 (x) 



2M~ l % m x 12 



-V5m + 2V3m) M p 2 x 10 



f A 2 „ \ o (2 4 \ ,, 2 fi 2 

+ ( g%?72 - g»?5»? - 2.gi?74 I x + I —77377 + -5i?72 I M pl x° + -giM pl r)x 



32g 4 - 12g 5 .3252 + 1253 4 
Mf, X + 



G(x) = (-^ +r?2X 4 + 7?M 2 x2 . ^ M li 



P i 



+ 30^A) 
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